The dynamics of non-linear optical absorption 
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On traversing materials with absorbing dyes, weak optical beams develop a Beer (exponential) 
profile, while intense beams develop a spatially initially linear and then finally an exponential profile. 
This anomalous, deep penetration due to photo-bleaching of surface layers is important for heavy 
dye-loading and intense beams, for instance in photo-actuation. We address the problem of the 
evolution in time from initial Beer's Law to deeply penetrating optical profiles in dyes. Our solution 
of the coupled, non-linear, partial differential equations governing the spatio-temporal decay of 
the Poynting flux and the non-linear dynamics of the trans-cis conversion is applicable to general 
systems of photo-active molecules under intense irradiation, for instance in biology, in spectroscopy 
and in opto-mechanical devices. 



Light absorption can optically switch the nematic 
phase of liquid crystals to the disordered, isotropic state: 
rodlike dye molecules transform from their linear trans 
to their bent cis isomers and thereby disrupt the orien- 
tational order of their nematic hosts. Although an inter- 
esting optical effect in its own right, when such a host is 
part of a nematic elastomer, then the solid responds with 
stresses Q or huge (100s%) optically-induced strains 
@; • Optical actuation offers advantages over electrical, 
thermal and solvent mediated response, being of large 
amplitude, easily reversible, sensitive to polarisation and 
susceptible to remote application. 

Two mysteries attend photoisomerisation, seen for in- 
stance in actuation. For linear absorptive processes at 
low light intensities (and therefore exponential profiles 
with depth, i.e. Beer's Law), and for mechanical contrac- 
tion proportional to trans-cis conversion, one can show 
4] that optimal bend of a photo-cantilever obtains at 
w/d ~ 3 where w is the cantilever thickness and d the 
penetration depth in Beer's Law 2 — I(x)/I(0) — e~ x / d 
with I(x) the Poynting flux at depth x in the cantilever. 
For iCw only a thin skin of material near the front face 
contracts and the remaining material resists bend. This 
limit of little penetration is that of high dye loading and 
in practice is often employed [Bj]. Non- linear absorption 
at high light intensities, where Beer's law gives way to a 
linear and more penetratingprofile for T(x) [a| , has been 
suggested as a mechanism p], @| to explain the bending 
of cantilevers with d -C w. The second mystery concerns 
the dynamical response, which is often initially slow and 
then proceeds quickly [9(] , especially in heavily dye-loaded 
systems. It is this second mystery that we address in this 
paper. 

In considering trans and cis populations and their dy- 
namics, we assume for simplicity that transitions are only 
optically-induced in the forward, t— >c direction and that 
the c— >t back-reaction is not optically, but only ther- 
mally induced (at rate 1/r), as is often the case [10(; see 
also [6( . We return to the case of an optically stimulated 



back-reaction briefly at the end of the paper. Given these 
optical transitions, the spatial decay of the Poynting flux 
of a light beam, I(x), is governed by 

dI/dx = -^Tnt{I)I=-n t I/d (1) 

where n t is the local number fraction of trans dye 
molecules and depends on time t and space x (through 
I). The Beer length is d — l/(-yT) that depends on the 
material parameters 7 (proportional to the number den- 
sity of chromophorcs and the energy each t— >c transition 
absorbs from the beam) and on T which determines the 
t— >c transition rate as we also see in the dynamical equa- 
tion for nt- 

dnt/dt = -TIrh + n c /r. (2) 

The cis number fraction is n c = 1 — nt- We reduce inten- 
sity by the incident value to give T(x) — I(x)/Iq. The 
combination I t = 1/(Tt) gives a characteristic intensity, 
a material constant related to the trans photo response. 
If Jo is reduced by It , then a = 1$ /It is a measure of how 
intense the incident beam is. The above equations then 
reduce to: 

dl/dx = -n t J/d , dnt/dt = -[{! + al)n t - 1] /t (3) 

At short times conversion has not yet proceeded and n t — 
1. Then the first of ([3|) easily integrates to Beer's Law 
l(x,t ~ 0) = e~ x l d . 

At long times (equilibrium) dnt /dt = and thus = 
1/(1 + olT). For weak beams (a ~ 0) n t « 1 and again 
([3]) gives Beer's Law. At high intensities a> 1 one has 
non-linear absorption [f|. Then n t ~ 1/(<xT), at least 
before depths such that the beam attenuates to I ~ I /a. 
Then in (J3J) one has dl/dx « — I /(ad) and hence X ~ 
1 — x/(ad) (for x < ad). The intensity profile is initially 
linear until depths where nt rises to being significantly 
greater than 1/a and absorption is then important. The 
surface layer is in effect photo-bleached and lets much 
light down to depths x > d; see @, H| for a detailed 
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discussion, including of the role of optical backreaction. 
We consider here for simplicity non mesogenic dyes and 
thus avoid dye rotation rather than cis formation as a 
route to bleaching (see ll| for this possibility) . With the 
equilibrium form of nt inserted into the first part of l[3"]). 
we can integrate to give the full, non-linear, equilibrium 
profile: 



ln[J(cc)] + a[I(x) - 1] = -x/d 



(4) 



The long-time limit (t > 5r in practice) of our dynami- 
cal profiles will display the linear rather than exponential 
forms until depths much greater than d. The Mathemat- 
ica defined function ProductLog gives 0, 0] the solution: 
1 = ^ProductLog [ae'"" 1 '] . 

We are concerned here with the dynamics of the tran- 
sition from Beer's law initially to the above equilib- 
rium profile for intense beams. The absorbance A = 
In [1/X(w, t)] is the usual measured quantity and reflects 
the absorber number in the optical path. Rearranging 
the first of eqns © to = —n t (x,t)/d and inte- 

grating J^dX and J™ dx, one obtains for all times and 
incident intensities (recall 2(0, t) — 1 for all t): 



A(w,t) = — n t 



w 1 
d w 



dx nt(x, t) 



(5) 



where fit is the mean trans number fraction through the 
sample. In the non- linear limit, fit is not independent 
of w and hence h t /d is no longer a simple, material- 
dependent extinction coefficient. 

The dynamics of A is often observed and analysed as- 
suming exponential behaviour: 

A(w, t) ~ A(w, oo) + [A(w,0) - A{w,oo)]e~ t/{n ~ T} . 

Note that A(w, 0) = w/d exactly and that the non-linear 
limit of A(w, oo) is given by the ProductLog solution to 
eqn ([4} for T(w) . However, this is not a solution to the 
differential eqns ((3|) and becomes a bad dynamical esti- 
mate for thick samples, w 3> d, that are only traversed 
by intense beams because of bleaching. The character- 
istic time tti%° = t/(1 + al) is shorter than thermal 
times when the driving a is large (intense beams), but 
clearly I(x) depends on position in thick samples and 
there is a spectrum of times and the overall response 
is not exponential. To solve the problem exactly, use 
A(x, t) = — ln(T) as the variable. Denote partial spatial 
and temporal derivatives of X by X' and X respectively. 
Differentiating eqn ([5|) with respect to t and using the 
second of eqns (j3|) for n t under the integral yields: 

A= [ dx [1 - (1 + al)n t ] /(rd) 
Jo 

Now use nt/d = A', A(0, t) = 0, and A 1 = -1' /I: 

tA = x/d- [ dx(l + al)A' (6) 
Jo 

= - ■- A + a(I - 1) = - ■- A + a (e~ A - l) .(7) 
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FIG. 1: Intensity against reduced depth for reduced incident 
intensities a — 10 and 30 (reduced times marked). 



A final quadrature gives A(x,t): 

dA 



t/r 



x/d 



x/d — a — A + ae 



-A' 



(8) 



The initial absorption A(x,t — 0) = x/d obtains from 
the vanishing of each side of ([8]) , and is Beer's Law. 

Figure [T] first shows X(x, t) (= e - " 4 ^'*)) for reduced in- 
tensity a = Iq I I c = 10 as a function of x for a set of times 
t. Initially at t = 0, the profile T(x, 0) is exponential, and 
at long times (t = 5r) the profile is essentially linear out 
to x ~ ad — lOd, then decays exponentially. At interme- 
diate times the profile first saturates (bleaches) at small 
x, that is it approaches the initial part of the equilibrium 
profile. Then as the surface layers let more light through, 
the profile deeper down also approaches the non-Beer 
form. For even higher a = 30, higher incident inten- 
sity, ultimate penetration is deeper and the approach to 
the bleached state, as time advances, even sharper. A 
front of bleaching propagates through the sample. Com- 
pare the two graphs at a given x — 5d, say. The more 
intense case, a = 30, sees a much quicker achievement of 
the bleached state than the a — 10 case, and over times 
considerably shorter than r, as the approximate analysis 
suggested. For either a, a sample with thickness w = 5d 7 
say, would see an intensity of emerging light far in excess 
of any Beer expectation, the rate build up of which we 
now examine. 
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FIG. 2: Emergent reduced intensity X = I(w)/Io against time 
at reduced incident intensity a = Io/Ic = 10 and 30 for vari- 
ous reduced sample thicknesses w/d. 



Although these curves of Fig.Q]are vital to understand- 
ing the build up (and possible decay) of optically-induced 
curvature in heavily dye-doped solid nematics, the inten- 
sity profile is not in general directly observable. It is 
more feasible to measure the dynamics of the build-up of 
intensity of the light emerging at the back face x = w of a 
sample, X(w,t). Figure reveals I(w, t) for various fixed 
thicknesses w, for two reduced incident intensities a. The 
initial (t = 0) value is just that from Beer penetration, 
e -w/d_ Finally T rises to the long-time, bleached value 
2(w,t ^> r) shown at the corresponding x = w in Figs[T] 
The rise is naturally slower for thicker samples. Such dy- 
namics of light penetration has been seen by Serra and 
Terentjev [9(. 

Thick (w > d) beams give initially exponentially small 
emergent fluxes that give way to large fluxes as the lin- 
ear profile is set up. Such data are accommodated in the 
usual logarithmic way, here by dividing the absorbance 
by the thickness which is not a material-dependent con- 
stant (an extinction) in this non- linear limit. Figure [3] 
shows A(w , t) I (w I d) for various fixed thicknesses w, for 
incident intensities a = 10, 30. 

The initial increase of intensity away from the Beer 
value in Fig. EU is j(w,t = 0) = -A{w,0)e' A ^ w '°'> = 
—A(w, 0)e~ w / d whereupon using for A at t = where 
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FIG. 3: Non-linear extinction A(w,t)/ (w/d) against time at 
reduced incident intensity a — Io/Ic = 10 and 30 for vari- 
ous reduced sample thicknesses w/d. The dotted line is the 
asymptote for the w/d = 30 dynamics. 

A = w/d, one obtains for initial slope 

T±{w,t = 0)=a(l-e- w/d ^je~ w/d . (9) 

For intense beams (a ^> 1) on thick (w ~ ad) samples, 
the initial intensities at the back face on irradiating the 
front face are small, X(w,0) — e~ w ^ d . The rates of in- 
crease are also initially very small, I ~ ae _Q /r, but then 
rise sharply with time, see the example with a = 30. 

Two other processes are sometimes important in 
non-linear absorption, those of host absorption and of 
cis absorption, that is optically-stimulated back reaction. 
Both are discussed in detail in the non-linear static case 
[8| . Moderate host absorption can be successfully divided 
out and plays a limited role. Back reaction occurs if the 
t— >c and c— >t absorption lines begin to overlap. This is 
sometimes the case and has the effect of reducing deep 
penetration 0]. 

In eqn ([1]) host absorption gives another term, —I/dh, 
while cis depletion of the beam gives — 7 c r c rt c (/)7. Op- 
tical back reaction also effects the dynamics. An addi- 
tional term +T c In c from the decay c— >t acts in eqn ([2]) 
to replenish the trans population. There is a charac- 
teristic intensity I c = l/(r c r) analogous to that of the 
trans species, and thus another measure (3 = Iq/I c of the 
incident beam intensity. Now the coupled, non-linear 
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FIG. 4: Emergent reduced light intensity T(w, t) against time 
for various reduced sample thicknesses w jd. Reduced incident 
intensity is a = Io/I c = 10 with a relatively high degree of 
reduced optical back reaction /3 = 2. 



partial differential equations (J3J) become: 

, _ dl/dx _ t\ _ l\ n l_ 
2 \d d c J d c 
rut = (1 + /3T) - [l + (a + /3)2]n t . 



1 

dh 



(10) 
(11) 



The changes, even if only one of these influences is intro- 
duced at a time, add difficulty to the solution of the equa- 
tions. Proceeding as before, one differentiates eqn (JTUJ) 
with respect to time and uses (fTTj) to eliminate n t . Where 
n t appears, it can be eliminated in favour of A' by using 
eqn (jTOj) again. Terms do not quite all integrate totally as 
in going from © to ([7]) because eqn (JTUJ) has additional 
new constant terms. Spatially integrating gives: 



tA= - A + x/d cB + (a + /3)(I-l) + 
1 1 
. d c dh 



dxl (12) 



where 1 / d cB = 1 / d+ 1 /dh is an effective absorption length 
arising from the simple combination of dye and host ab- 
sorptions in the Beer limit. It is interesting that although 
f} (and thus l/d c ) and l/dh enter eqns (fTOj) and (fTTj) in 
entirely different ways, their effects can be scaled on to 
each other in the resultant equation for dynamical non- 
linear absorption (|12[) . Remember however that d e ff de- 
pends on dh and so lengths are not quite equivalcntly 



effected by this interchange of sources of extra absorp- 
tion. We show just the effect of adding in optical back- 
reaction into the dynamical equations, that is 1/d^ = 
and lengths are still reduced by d. For a = 10 and, say, 
(3 = 0.2, the I{w, t)—t curves are indistinguishable from 
those in upper Fig. [5] which have (3 — 0. Figure Q] shows 
the behaviour for f3 — 2 which should be compared with 
upper Fig. [2] The starting values are of course the same; 
the final values are lower because j3 acts to lower the final 
penetration. The upward curvature at short times is lost. 
Thicker samples are much more drastically affected. 

We have shown that the dynamics of penetration of 
intense light beams into heavily absorbing media is com- 
plex and strongly dependent upon incident intensity. A 
practical experiment to observe this effect via the emer- 
gent light from the sample is proposed. There is evidence 
from large amplitude photo-mechanics that non-linear ef- 
fects are important in practice. 
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